Single-cell and spatial transcriptomics identify a macrophage population associated with skeletal muscle fibrosis.

Macrophages are essential for skeletal muscle homeostasis, but how their dysregulation contributes to the development of fibrosis in muscle disease remains unclear. Here, we used single-cell transcriptomics to determine the molecular attributes of dystrophic and healthy muscle macrophages. We identified six clusters and unexpectedly found that none corresponded to traditional definitions of M1 or M2 macrophages. Rather, the predominant macrophage signature in dystrophic muscle was characterized by high expression of fibrotic factors, galectin-3 (gal-3) and osteopontin (Spp1). Spatial transcriptomics, computational inferences of intercellular communication, and in vitro assays indicated that macrophage-derived Spp1 regulates stromal progenitor differentiation. Gal-3+ macrophages were chronically activated in dystrophic muscle, and adoptive transfer assays showed that the gal-3+ phenotype was the dominant molecular program induced within the dystrophic milieu. Gal-3+ macrophages were also elevated in multiple human myopathies. These studies advance our understanding of macrophages in muscular dystrophy by defining their transcriptional programs and reveal Spp1 as a major regulator of macrophage and stromal progenitor interactions.

Single-cell and spatial transcriptomics identify a macrophage population associated with skeletal muscle fibrosis.

The monocytic/macrophage system is essential for skeletal muscle homeostasis, but its dysregulation contributes to the pathogenesis of muscle degenerative disorders. Despite our increasing knowledge of the role of macrophages in degenerative disease, it still remains unclear how macrophages contribute to muscle fibrosis. Here, we used single-cell transcriptomics to determine the molecular attributes of dystrophic and healthy muscle macrophages. We identified six novel clusters. Unexpectedly, none corresponded to traditional definitions of M1 or M2 macrophage activation. Rather, the predominant macrophage signature in dystrophic muscle was characterized by high expression of fibrotic factors, galectin-3 and spp1. Spatial transcriptomics and computational inferences of intercellular communication indicated that spp1 regulates stromal progenitor and macrophage interactions during muscular dystrophy. Galectin-3 + macrophages were chronically activated in dystrophic muscle and adoptive transfer assays showed that the galectin-3 + phenotype was the dominant molecular program induced within the dystrophic milieu. Histological examination of human muscle biopsies revealed that galectin-3 + macrophages were also elevated in multiple myopathies. These studies advance our understanding of macrophages in muscular dystrophy by defining the transcriptional programs induced in muscle macrophages, and reveal spp1 as a major regulator of macrophage and stromal progenitor interactions.

Overexpressed Malat1 Drives Metastasis through Inflammatory Reprogramming of Lung Adenocarcinoma Microenvironment.

Metastasis is the main cause of cancer deaths but the molecular events leading to metastatic dissemination remain incompletely understood. Despite reports linking aberrant expression of long noncoding RNAs (lncRNAs) with increased metastatic incidence , in vivo evidence establishing driver roles for lncRNAs in metastatic progression is lacking. Here, we report that overexpression of the metastasis-associated lncRNA Malat1 (metastasis-associated lung adenocarcinoma transcript 1) in the autochthonous K-ras/p53 mouse model of lung adenocarcinoma (LUAD) is sufficient to drive cancer progression and metastatic dissemination. We show that increased expression of endogenous Malat1 RNA cooperates with p53 loss to promote widespread LUAD progression to a poorly differentiated, invasive, and metastatic disease. Mechanistically, we observe that Malat1 overexpression leads to the inappropriate transcription and paracrine secretion of the inflammatory cytokine, Ccl2, to augment the mobility of tumor and stromal cells in vitro and to trigger inflammatory responses in the tumor microenvironment in vivo . Notably, Ccl2 blockade fully reverses cellular and organismal phenotypes of Malat1 overexpression. We propose that Malat1 overexpression in advanced tumors activates Ccl2 signaling to reprogram the tumor microenvironment to an inflammatory and pro-metastatic state.

The p53 transcriptional response across tumor types reveals core and senescence-specific signatures modulated by long noncoding RNAs.

The p53 pathway is a universal tumor suppressor mechanism that limits tumor progression by triggering apoptosis or permanent cell cycle arrest, called senescence. In recent years, efforts to reactivate p53 function in cancer have proven to be a successful therapeutic strategy in murine models and have gained traction with the development of a range of small molecules targeting mutant p53. However, knowledge of the downstream mediators of p53 reactivation in different oncogenic contexts has been limited. Here, we utilized a panel of murine cancer cell lines from three distinct tumor types susceptible to alternative outcomes following p53 restoration to define unique and shared p53 transcriptional signatures. While we found that the majority of p53-bound sites and p53-responsive transcripts are tumor-type specific, analysis of shared targets identified a core signature of genes activated by p53 across all contexts. Furthermore, we identified repression of E2F and Myc target genes as a key feature of senescence. Characterization of p53-induced transcripts revealed core and senescence-specific long noncoding RNAs (lncRNAs) that are predominantly chromatin associated and whose production is coupled to cis-regulatory activities. Functional investigation of the contributions of p53-induced lncRNAs to p53-dependent outcomes highlighted Pvt1b, the p53-dependent isoform of Pvt1, as a mediator of p53-dependent senescence via Myc repression. Inhibition of Pvt1b led to decreased activation of senescence markers and increased levels of markers of proliferation. These findings shed light on the core and outcome-specific p53 restoration signatures across different oncogenic contexts and underscore the key role of the p53-Pvt1b-Myc regulatory axis in mediating proliferative arrest.

p53 Activates the Long Noncoding RNA Pvt1b to Inhibit Myc and Suppress Tumorigenesis.

The tumor suppressor p53 transcriptionally activates target genes to suppress cellular proliferation during stress. p53 has also been implicated in the repression of the proto-oncogene Myc, but the mechanism has remained unclear. Here, we identify Pvt1b, a p53-dependent isoform of the long noncoding RNA (lncRNA) Pvt1, expressed 50 kb downstream of Myc, which becomes induced by DNA damage or oncogenic signaling and accumulates near its site of transcription. We show that production of the Pvt1b RNA is necessary and sufficient to suppress Myc transcription in cis without altering the chromatin organization of the locus. Inhibition of Pvt1b increases Myc levels and transcriptional activity and promotes cellular proliferation. Furthermore, Pvt1b loss accelerates tumor growth, but not tumor progression, in an autochthonous mouse model of lung cancer. These findings demonstrate that Pvt1b acts at the intersection of the p53 and Myc transcriptional networks to reinforce the anti-proliferative activities of p53.

Temporal Control of the TGF-ß Signaling Network by Mouse ESC MicroRNA Targets of Different Affinities.

Although microRNAs (miRNAs) function in the control of embryonic stem cell (ESC) pluripotency, a systems-level understanding is still being developed. Through the analysis of progressive Argonaute (Ago)-miRNA depletion and rescue, including stable Ago knockout mouse ESCs, we uncover transforming growth factor beta (TGF-ß) pathway activation as a direct and early response to ESC miRNA reduction. Mechanistically, we link the derepression of weaker miRNA targets, including TGF-ß receptor 1 (Tgfbr1), to the sensitive TGF-ß pathway activation. In contrast, stronger miRNA targets impart a more robust repression, which dampens concurrent transcriptional activation. We verify such dampened induction for TGF-ß antagonist Lefty. We find that TGF-ß pathway activation contributes to the G1 cell-cycle accumulation of miRNA-deficient ESCs. We propose that miRNA target affinity is a determinant of the temporal response to miRNA changes, which enables the coordination of gene network responses.

Sequestration of microRNA-mediated target repression by the Ago2-associated RNA-binding protein FAM120A.

Argonaute (Ago) proteins interact with various binding partners and play a pivotal role in microRNA (miRNA)-mediated silencing pathways. By utilizing immunoprecipitation followed by mass spectrometry to determine cytoplasmic Ago2 protein complexes in mouse embryonic stem cells (mESCs), we identified a putative RNA-binding protein FAM120A (also known as OSSA/C9ORF10) as an Ago2 interacting protein. Individual nucleotide resolution cross-linking and immunoprecipitation (iCLIP) analysis revealed that FAM120A binds to homopolymeric tracts in 3'-UTRs of about 2000 mRNAs, particularly poly(G) sequences. Comparison of FAM120A iCLIP and Ago2 iCLIP reveals that greater than one-third of mRNAs bound by Ago2 in mESCs are co-bound by FAM120A. Furthermore, such FAM120A-bound Ago2 target genes are not subject to Ago2-mediated target degradation. Reporter assays suggest that the 3'-UTRs of several FAM120A-bound miRNA target genes are less sensitive to Ago2-mediated target repression than those of FAM120A-unbound miRNA targets and FAM120A modulates them via its G-rich target sites. These findings suggest that Ago2 may exist in multiple protein complexes with varying degrees of functionality.

Endogenous miRNA and target concentrations determine susceptibility to potential ceRNA competition.

Target competition (ceRNA crosstalk) within miRNA-regulated gene networks has been proposed to influence biological systems. To assess target competition, we characterize and quantitate miRNA networks in two cell types. Argonaute iCLIP reveals that hierarchical binding of high- to low-affinity miRNA targets is a key characteristic of in vivo activity. Quantification of cellular miRNA and mRNA/ncRNA target pool levels indicates that miRNA:target pool ratios and an affinity partitioned target pool accurately predict in vivo Ago binding profiles and miRNA susceptibility to target competition. Using single-cell reporters, we directly test predictions and estimate that ?3,000 additional high-affinity target sites can affect active miRNA families with low endogenous miRNA:target ratios, such as miR-92/25. In contrast, the highly expressed miR-294 and let-7 families are not susceptible to increases of nearly 10,000 sites. These results show differential susceptibility based on endogenous miRNA:target pool ratios and provide a physiological context for ceRNA competition in vivo.

LincRNA-p21 activates p21 in cis to promote Polycomb target gene expression and to enforce the G1/S checkpoint.

The p53-regulated long noncoding RNA lincRNA-p21 has been proposed to act in trans via several mechanisms ranging from repressing genes in the p53 transcriptional network to regulating mRNA translation and protein stability. To further examine lincRNA-p21 function, we generated a conditional knockout mouse model. We find that lincRNA-p21 predominantly functions in cis to activate expression of its neighboring gene, p21. Mechanistically, we show that lincRNA-p21 acts in concert with hnRNP-K as a coactivator for p53-dependent p21 transcription. Additional phenotypes of lincRNA-p21 deficiency could be attributed to diminished p21 levels, including deregulated expression and altered chromatin state of some Polycomb target genes, a defective G1/S checkpoint, increased proliferation rates, and enhanced reprogramming efficiency. These findings indicate that lincRNA-p21 affects global gene expression and influences the p53 tumor suppressor pathway by acting in cis as a locus-restricted coactivator for p53-mediated p21 expression.

Argonaute-bound small RNAs from promoter-proximal RNA polymerase II.

Argonaute (Ago) proteins mediate posttranscriptional gene repression by binding guide miRNAs to regulate targeted RNAs. To confidently assess Ago-bound small RNAs, we adapted a mouse embryonic stem cell system to express a single epitope-tagged Ago protein family member in an inducible manner. Here, we report the small RNA profile of Ago-deficient cells and show that Ago-dependent stability is a common feature of mammalian miRNAs. Using this criteria and immunopurification, we identified an Ago-dependent class of noncanonical miRNAs derived from protein-coding gene promoters, which we name transcriptional start site miRNAs (TSS-miRNAs). A subset of promoter-proximal RNA polymerase II (RNAPII) complexes produces hairpin RNAs that are processed in a DiGeorge syndrome critical region gene 8 (Dgcr8)/Drosha-independent but Dicer-dependent manner. TSS-miRNA activity is detectable from endogenous levels and following overexpression of mRNA constructs. Finally, we present evidence of differential expression and conservation in humans, suggesting important roles in gene regulation.

Antisense RNA polymerase II divergent transcripts are P-TEFb dependent and substrates for the RNA exosome.

Divergent transcription occurs at the majority of RNA polymerase II (RNAPII) promoters in mouse embryonic stem cells (mESCs), and this activity correlates with CpG islands. Here we report the characterization of upstream antisense transcription in regions encoding transcription start site associated RNAs (TSSa-RNAs) at four divergent CpG island promoters: Isg20l1, Tcea1, Txn1, and Sf3b1. We find that upstream antisense RNAs (uaRNAs) have distinct capped 5' termini and heterogeneous nonpolyadenylated 3' ends. uaRNAs are short-lived with average half-lives of 18 minutes and are present at 1-4 copies per cell, approximately one RNA per DNA template. Exosome depletion stabilizes uaRNAs. These uaRNAs are probably initiation products because their capped termini correlate with peaks of paused RNAPII. The pausing factors NELF and DSIF are associated with these antisense polymerases and their sense partners. Knockdown of either NELF or DSIF results in an increase in the levels of uaRNAs. Consistent with P-TEFb controlling release from pausing, treatment with its inhibitor, flavopiridol, decreases uaRNA and nascent mRNA transcripts with similar kinetics. Finally, Isg20l1 induction reveals equivalent increases in transcriptional activity in sense and antisense directions. Together these data show divergent polymerases are regulated after P-TEFb recruitment with uaRNA levels controlled by the exosome.

2'-O-ribose methylation of cap2 in human: function and evolution in a horizontally mobile family.

The 5' cap of human messenger RNA consists of an inverted 7-methylguanosine linked to the first transcribed nucleotide by a unique 5'-5' triphosphate bond followed by 2'-O-ribose methylation of the first and often the second transcribed nucleotides, likely serving to modify efficiency of transcript processing, translation and stability. We report the validation of a human enzyme that methylates the ribose of the second transcribed nucleotide encoded by FTSJD1, henceforth renamed HMTR2 to reflect function. Purified recombinant hMTr2 protein transfers a methyl group from S-adenosylmethionine to the 2'-O-ribose of the second nucleotide of messenger RNA and small nuclear RNA. Neither N(7) methylation of the guanosine cap nor 2'-O-ribose methylation of the first transcribed nucleotide are required for hMTr2, but the presence of cap1 methylation increases hMTr2 activity. The hMTr2 protein is distributed throughout the nucleus and cytosol, in contrast to the nuclear hMTr1. The details of how and why specific transcripts undergo modification with these ribose methylations remains to be elucidated. The 2'-O-ribose RNA cap methyltransferases are present in varying combinations in most eukaryotic and many viral genomes. With the capping enzymes in hand their biological purpose can be ascertained.

Hypermethylated cap 4 maximizes Trypanosoma brucei translation.

Through trans-splicing of a 39-nt spliced leader (SL) onto each protein-coding transcript, mature kinetoplastid mRNA acquire a hypermethylated 5'-cap structure, but its function has been unclear. Gene deletions for three Trypanosoma brucei cap 2'-O-ribose methyltransferases, TbMTr1, TbMTr2 and TbMTr3, reveal distinct roles for four 2'-O-methylated nucleotides. Elimination of individual gene pairs yields viable cells; however, attempts at double knock-outs resulted in the generation of a TbMTr2-/-/TbMTr3-/- cell line only. Absence of both kinetoplastid-specific enzymes in TbMTr2-/-/TbMTr3-/- lines yielded substrate SL RNA and mRNA with cap 1. TbMTr1-/- translation is comparable with wildtype, while cap 3 and cap 4 loss reduced translation rates, exacerbated by the additional loss of cap 2. TbMTr1-/- and TbMTr2-/-/TbMTr3-/- lines grow to lower densities under normal culture conditions relative to wildtype cells, with growth rate differences apparent under low serum conditions. Cell viability may not tolerate delays at both the nucleolar Sm-independent and nucleoplasmic Sm-dependent stages of SL RNA maturation combined with reduced rates of translation. A minimal level of mRNA cap ribose methylation is essential for trypanosome viability, providing the first functional role for the cap 4.

Trypanosoma brucei spliced leader RNA maturation by the cap 1 2'-O-ribose methyltransferase and SLA1 H/ACA snoRNA pseudouridine synthase complex.

Kinetoplastid flagellates attach a 39-nucleotide spliced leader (SL) upstream of protein-coding regions in polycistronic RNA precursors through trans splicing. SL modifications include cap 2'-O-ribose methylation of the first four nucleotides and pseudouridine (psi) formation at uracil 28. In Trypanosoma brucei, TbMTr1 performs 2'-O-ribose methylation of the first transcribed nucleotide, or cap 1. We report the characterization of an SL RNA processing complex with TbMTr1 and the SLA1 H/ACA small nucleolar ribonucleoprotein (snoRNP) particle that guides SL psi(28) formation. TbMTr1 is in a high-molecular-weight complex containing the four conserved core proteins of H/ACA snoRNPs, a kinetoplastid-specific protein designated methyltransferase-associated protein (TbMTAP), and the SLA1 snoRNA. TbMTAP-null lines are viable but have decreased SL RNA processing efficiency in cap methylation, 3'-end maturation, and psi(28) formation. TbMTAP is required for association between TbMTr1 and the SLA1 snoRNP but does not affect U1 small nuclear RNA methylation. A complex methylation profile in the mRNA population of TbMTAP-null lines indicates an additional effect on cap 4 methylations. The TbMTr1 complex specializes the SLA1 H/ACA snoRNP for efficient processing of multiple modifications on the SL RNA substrate.

The TbMTr1 spliced leader RNA cap 1 2'-O-ribose methyltransferase from Trypanosoma brucei acts with substrate specificity.

In metazoa cap 1 (m(7)GpppNmp-RNA) is linked to higher levels of translation; however, the enzyme responsible remains unidentified. We have validated the first eukaryotic encoded cap 1 2'-O-ribose methyltransferase, TbMTr1, a member of a conserved family that modifies the first transcribed nucleotide of spliced leader and U1 small nuclear RNAs in the kinetoplastid protozoan Trypanosoma brucei. In addition to cap 0 (m(7)GpppNp-RNA), mRNA in these parasites has ribose methylations on the first four nucleotides with base methylations on the first and fourth (m(7)Gpppm(6,6)AmpAmpCmpm(3)Ump-SL RNA) conveyed via trans-splicing of a universal spliced leader. The function of this cap 4 is unclear. Spliced leader is the majority RNA polymerase II transcript; the RNA polymerase III-transcribed U1 small nuclear RNA has the same first four nucleotides as spliced leader, but it receives an m(2,2,7)G cap with hypermethylation at position one only (m(2,2,7)Gpppm(6,6)AmpApCpUp-U1 snRNA). Here we examine the biochemical properties of recombinant TbMTr1. Active over a pH range of 6.0 to 9.5, TbMTr1 is sensitive to Mg(2+). Positions Lys(95)-Asp(204)-Lys(259)-Glu(285) constitute the conserved catalytic core. A guanosine cap on RNA independent of its N(7) methylation status is required for substrate recognition, but an m(2,2,7G)-cap is not recognized. TbMTr1 favors the spliced leader 5' sequence, as reflected by a preference for A at position 1 and modulation of activity for substrates with base changes at positions 2 and 3. With similarities to human cap 1 methyltransferase activity, TbMTr1 is an excellent model for higher eukaryotic cap 1 methyltransferases and the consequences of cap 1 modification.

The 2'-O-ribose methyltransferase for cap 1 of spliced leader RNA and U1 small nuclear RNA in Trypanosoma brucei.

mRNA cap 1 2'-O-ribose methylation is a widespread modification that is implicated in processing, trafficking, and translational control in eukaryotic systems. The eukaryotic enzyme has yet to be identified. In kinetoplastid flagellates trans-splicing of spliced leader (SL) to polycistronic precursors conveys a hypermethylated cap 4, including a cap 0 m7G and seven additional methylations on the first 4 nucleotides, to all nuclear mRNAs. We report the first eukaryotic cap 1 2'-O-ribose methyltransferase, TbMTr1, a member of a conserved family of viral and eukaryotic enzymes. Recombinant TbMTr1 methylates the ribose of the first nucleotide of an m7G-capped substrate. Knockdowns and null mutants of TbMTr1 in Trypanosoma brucei grow normally, with loss of 2'-O-ribose methylation at cap 1 on substrate SL RNA and U1 small nuclear RNA. TbMTr1-null cells have an accumulation of cap 0 substrate without further methylation, while spliced mRNA is modified efficiently at position 4 in the absence of 2'-O-ribose methylation at position 1; downstream cap 4 methylations are independent of cap 1. Based on TbMTr1-green fluorescent protein localization, 2'-O-ribose methylation at position 1 occurs in the nucleus. Accumulation of 3'-extended SL RNA substrate indicates a delay in processing and suggests a synergistic role for cap 1 in maturation.

Complete cap 4 formation is not required for viability in Trypanosoma brucei.

In kinetoplastids spliced leader (SL) RNA is trans-spliced onto the 5' ends of all nuclear mRNAs, providing a universal exon with a unique cap. Mature SL contains an m(7)G cap, ribose 2'-O methylations on the first four nucleotides, and base methylations on nucleotides 1 and 4 (AACU). This structure is referred to as cap 4. Mutagenized SL RNAs that exhibit reduced cap 4 are trans-spliced, but these mRNAs do not associate with polysomes, suggesting a direct role in translation for cap 4, the primary SL sequence, or both. To separate SL RNA sequence alterations from cap 4 maturation, we have examined two ribose 2'-O-methyltransferases in Trypanosoma brucei. Both enzymes fall into the Rossmann fold class of methyltransferases and model into a conserved structure based on vaccinia virus homolog VP39. Knockdown of the methyltransferases individually or in combination did not affect growth rates and suggests a temporal placement in the cap 4 formation cascade: TbMT417 modifies A(2) and is not required for subsequent steps; TbMT511 methylates C(3), without which U(4) methylations are reduced. Incomplete cap 4 maturation was reflected in substrate SL and mRNA populations. Recombinant methyltransferases bind to a methyl donor and show preference for m(7)G-capped RNAs in vitro. Both enzymes reside in the nucleoplasm. Based on the cap phenotype of substrate SL stranded in the cytosol, A(2), C(3), and U(4) methylations are added after nuclear reimport of Sm protein-complexed substrate SL RNA. As mature cap 4 is dispensable for translation, cap 1 modifications and/or SL sequences are implicated in ribosomal interaction.